Error Bounds for Asymptotic Solutions of 
Differential Equations. 
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The results of the preceding paper are extended to the general system of n first-order differential 
equations having an irregular singularity of arbitrary rank at infinity. Formal solutions are explicitly 
constructed for the system in canonical form. Proofs of existence and uniqueness of solutions of 
integral equations defining the error are given. As an example, the case n = 2 is solved completely, 
and a flow chart of the transformations of this case to canonical form is included. 



1. Introduction 

The system of differential equations under consideration in this paper is of the form 

W, 5 = 0, 1, . . ., (1.1) 



dz 



A-=0 



where h is an integer and the matrices A* and ^ls(z) are defined as in [1], except that no assump- 
tions are made concerning eigenvalues of the matrices A* (& = 0, 1, . . .)• We shall simulta- 
neously consider the case of a regular point, the case of a regular singular point, and the case of 
an irregular singular point at infinity. 

From the point of view of applications it is worthwhile recalling that when the eigenvalues of 
the lead coefficient matrix A () are not distinct they can be made distinct by an arbitrarily small 
change in the elements of A . The results of fl] 3 can then be applied, and the very complicated 
procedure given in this paper for obtaining a formal solution of (1.1) avoided. Indeed the ex- 
istence of a linearly independent set of n formal solution vectors for the system (1.1) was fully 
established only in 1937 ([2]) and more explicitly in 1955 ([3]). 

We have not been able to improve on the method of [3J for the reduction of (1.1) to canonical 
form. We thus start with the canonical form given in [3] and establish a minor variation of it 
which enables us more easily to obtain error bounds. Next we explicitly construct formal solu- 
tions for the system in canonical form. Integral equations are then set up for the error vector, 
that is, for the difference between the partial sum of a formal solution vector and an actual solution 
vector. 

As in [1], it is again possible to express the error by one or, at most, two Volterra vector integral 
equations. As opposed to [1], however, it is not always possible to choose the end points of in- 
tegration in the Volterra integral equations to be at infinity. 

2. The Canonical Form 

The following theorem due to H. L. Turrittin [3] concerns transformation to canonical form, of 
the system (1.1), with eigenvalues of A not necessarily distinct. In this theorem the normalizing, 
exponential, zero-inducing, root-equalizing, and shearing transformations are defined as in [3]. 
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An example of each of these transformations will be given in our treatment of the special n = 2 
case in section 4. 

THEOREM 2.1: Let h be an integer and 



dX 
dt ' 



t h ?l(t)X 



(2.1) 



a system of differential equations in which 

9(t)-§ A k r k + 2«t)t- s (s = 0, 1, . . .), 



(2.2) 



where the A k are constant n X n matrices and the elements ofilJX) (s = l,2, . . .) are holomorphic 
and uniformly bounded in some domain 3>\ extending to infinity. Then the system (2.1) can be 
transformed by a finite sequence of normalizing, zero-inducing, root-equalizing and shearing trans- 
formations to the canonical system 



where z = t p , r and p > are integers, 



dW 

^=z r G(z)W, 
dz 



(2.3) 



r (i(z) = 



e,(z) 



©2(2 



e,(z) 



+ z~ 2 



93„(z)93 12 (z) . . . ®«(z) 
93 2 i(z)93 22 (z) . . . ©„(z) 
«n(z)»«(a) • • • »«(z) 



(2.4) 



z r 6 1 (z)=2C ik z r - k I M . + z-'J M . 

k=0 

i=l,2, ...,/; /LtiH-/X2+. . .-\-/JLi=n 

«ij(z) = S f B ijk z- k + 23 ijs (z)z- s , 



(2.5) 



The elements of the matrices 33ij(z) are holomorphic for z in the domain 3) which corresponds to 
3) i for all z sufficiently large. Each element o/ 33 ijs (z) (i, j = l, 2, . . ., /; s = 0, 1, . . .) is uni- 
formly bounded on any path @> extending to x in 3 '. Ifr^—2 the matrices &i(z) (i = l, 2, . . . /) 
are identically zero. Ifr>—2, the l^.are /XiX/^ unit matrices, and the J M . are /jl { X ^ square matrices 
with zeros or ones or a mixture of zeros and ones in the first subdiagonal while all other elements of 
J^ are zero. No two of the matrices Gj(z) (i= 1, 2, . . ., I) in (2.5) are identical. In particular, 
if i # j and c ik — c ik for k = 0, 1 , . . ., r, then corresponding to any particular fixed] (j = 1, 2, . . . , /) 
we have — 1 ^ Rec h r+ i < 0, ^ ReCj, r +i < 1 and c u r +i — c L r+ i # 1 (i = l, 2, . . ., /; i # j). 

The statement of the above theorem is different in two aspects from that given by Turrittin 
in [3]. In the first place Turrittin proves the theorem for the case when ?((t) has a convergent 
power series expansion about t = 0. Secondly, Turrittin's version of the last sentence of the 
theorem is: In particular, i/i# j and c ik =c jk /or k = 0, 1, . . ., r then c u r+1 — c L T+1 is not only 
not zero but is also not an integer. 
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It is readily verified that all of the transformations: normalizing, zero-inducing, root-equalizing, 
and shearing preserve holomorphy of ?l(£) for all t sufficiently large, and £Ei^extending to infinity. 

We now establish the statement in the last sentence of Theorem 2.1, assuming Turrittin's 
version of the theorem. Consider the subset of integers {/} of the set {i = l, 2, . . ., /} such 
that with iG{i}, and i 9^ j\ c;a- — Cja = 0, k := 0, 1, . . ., r. Corresponding to the points Ci, r +i in 
the complex A plane, we define integers kj as follows: if a, r +i lies in the strip k% ^ Re \<k\+\ we 
associate the integer k\ with a, r +i. Without loss of generality we assume that {i} ={1, 2, . . ., a>} 
and that k\=* k 2 = . . . = A^. Let k = k lo — k\. 

We apply the sequence of zero-inducing transformations 



Z*+i = (I + TjbZ-*)Z*; A=l,2, . . .,/c+l; Z 2 = W 



(2.6) 



to (2.3), where TV is an n X n matrix with an arbitrary leading at X co block and zeros everywhere 
else. Thus for example if 



dZ k 



dz 



■=z r [Co+Ci2r 1 + . • . + C,. + iz- / - 1 + B () z- / - 



■ + B,,z~ 



[ + , 



.]Z A - 



(2.7) 



where the leading wXw blocks of B , Bi, . . ., Ba_! have already been filled with zeros, then 



dZ k _ 

dz 



■z r [C +Ciz- l + . . . + C r+1 2-'- | + B o2 -'- 2 + . . . + B*z- fc -'- 1 + . . .]Z* 



(2.8) 



yhere 



B* - B, - T*C r+1 + C r+ iT A . + /T,. 



(2.9) 



We divide the leading ojXoj blocks of the constant matrices B*, B&, Ta and C r +i into smaller sub- 
blocks corresponding to the way in which B () was subdivided in the above theorem. We denote 
these subblocks respectively by B?- A , By*, Ty&, and Cy, r +i- Equation (2.9) for Bf jk then becomes 



B 



W ~ 



By* + (C J M . H- J/A^Tyfe — T/ja-(CjI Mj + J Mj ) H- /cTjja-; 



(2.10) 



where C/ = C,-, r +i. Assuming Turrittin's version of Theorem 1, we have Xj— ki^ k. Elements 
of T/jA can therefore be computed 4 starting from the right-hand vector, calculating from top to 
bottom and in this manner successively computing the vectors from right to left such that B,* A =0. 
It is clear that each leading <o X ay block of B* (Ar = 0, 1, . . . , k) can be filled with zeros in this 
manner. 

The root-equalizing transformations 



z^L 



Z A ' 2 I 



Y, 



Y = 



1, 



**i-i 



zl 



*H 



z k l u 



are now applied to (2.8) above, changing Cj, r +i (t = l, 2, . 
sentence of Theorem 2.1. 



W 



/) to the form described in the last 



4 For further details, see [3] p. 33. 
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(2.12) 



The following theorem is also proved by Turrittin [3]. The main purpose of presenting a 
proof here is to obtain a recurrence relation for the coefficients of the solution matrix and to furnish 
an explicit algorithm for computing these coefficients. 

THEOREM 2.2: The system of differential equations in the canonical form (2.3) to (2.5) possesses 
a formal independent series solution matrix of the form 

W(z)=U(z)e c « (2.11) 

where 

U(z) = [U u (z)] ± f U k z- k * J [U iJk ]z" k 

k=0 k=0 

Uij(z) is /Xj X l± h Uyo= Syl^ 

G(z) = [(I^qrfz) + J^ log z)8 u ] = V <fr z r + ^ + Q r+ , log z 

k=o rtl R 

qi( z ) = i) r+ c i 11 Li f 2r+1 " k +qi,r + iiogz a,j=i,2, . . ., /). 

k=() r-ri k 
Proof of Theorem 2.2: Let us put 

z r [_Uz%j\ = 2 C^* 1 "* - z'D(z), (2.13) 

[»«&] = 23(z) = 2 B ft 2T* + »£)r-. (2.14) 

We substitute (2.11) into (2.3), and use the fact that Q and G' commute, to find 

U'(z) + U(z)C'(z) - [z>'D(z) 4- z- 2 93(z)]U(z). (2.15) 

Expanding (2.3), using (2.12), (2.13), and (2.14) and equating equal powers of z, we obtain 

min (u, r+1) fx-r-2 

2 (C fc U M _ fc -U M - fc Cit) + (/i-r-l)U M - r -,+ ^ B fc U M _ r _ 2 -* = (2.16) 

fc=0 Ar=0 



for /x = 0, 1, 2, . . ., where U&=Oif A<0 and C* = Q&=0 if *>r+l. 

The proof of Theorem 2.2 will be complete once we prove the following lemma. 
Lemma 2.1: Under the assumption that O(z) has the form given by (2.12) with Uo = I, equation 
(2.16) uniquely determines Q M (^ = 0, 1, 2, . . . r + 1) and U M (/n=l, 2, 3, . . .). 

Proof of Lemma 2.1: The proof of the lemma is by induction. With fi = 0, (2.16) gives 

C -Qo = 0. (2.17) 



which has the unique solution Qo = Co. With /x=l, (2.15) becomes 

Ci-Q 1 -hCoUi-UiC = 0. (2.18) 
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Now if Ui is an arbitrary n X n matrix we know from the definition of C that the (i,i)th (i = 1,2, 
...,/) blocks of Colli — U1C0 consisting of jjh X ^ matrices are identically zero. Hence these 
blocks of Ci — Qi must also be zero, and therefore, Qi = Ci. Furthermore, it follows that all (i,y)th 
blocks of Ui for which c io — Cjo ¥* are uniquely determined and are zero. All (i,y')th blocks of U, 
for which c i0 = cp cannot be determined at this point. In particular, it is not possible to solve for 

u«,o:=i.2. . . .,/). 

Assume that for some integer fe, ^ h = r, we have completed solving (2.16) for all /JL = k; 
that we have found Q M — C^ (/x = 0, 1, . . . , k) and that we have determined all blocks of the mat- 
rices U M as far as possible. More precisely, we make the following assumptions for each 

(i,/)th block of U M (i,y=l, 2, ...,/; /x = 0, 1 A). If c itl -c jtl = (/ll=0, 1, . . . , k-\) 

then assume U#i, U/j 2 , . . • , U#* are arbitrary, /JLiX/JLj matrices to be determined later. On the 
other hand, if /; is the smallest integer in ^ p ^ k for which Ctp—cjp 7^ 0, we assume that Utp 
= Uyi=. . . = U/j, a--/>— 0(/^/X/xj), while Uy, t+i-p, • • ■ , U#* are arbitrary /x/X^j matrices. 

For /jl = k-\- 1, equation (2.16) becomes 

CfU» + i-U* + iCo+. • . + C a -U 1 -U 1 C,+ C a+ ,-0a + i = O. (2.19) 

If Ui, U 2 , . . ., Ua + i are arbitrary a X n matrices, the (£, t)th (i= 1, 2, . . ., /) blocks of C^Ufc+i-M 
— U*+i- M C M are zero when /Lt=0, 1, . . . , A. Hence these blocks of Cfc+i — Qa+i must also be 
zero, and therefore Q*+i = C*+i. On incorporating this result in (2.19) we find that the (i,/)th 
block of (2.19) becomes 

(Cio — C/o)Uy,ifc + i + (Cji — Cji)Uy,*+. • .+(Cifc — C;*)U(ji = 0. (2.20) 

If c,> — Cj M = 0, /x=0, 1, 2, ... , k, then we cannot solve for U/ ; i, U#2, • • . , U/j, a + l We leave 
these matrices undetermined for the moment. On the other hand, if/; is the smallest nonnegative 
integer in ^ p ^ k for which dp — Cjp ^ 0, then by the above induction hypothesis, (2.20) becomes 
(cip — Cj P )\Jij t a-+i-/,=iO(/z,/X/zj), and so under the assumption that Uy,*+2-p, . . . , Uy,*+i are arbi- 
trary fjiiXfjij matrices we have Uy, a + i-/, = 0. 

It is clear at this point that Q M is uniquely determined for each /ll = 0, 1, . . . , r+ 1. We may 
therefore write the (i, /)th block of (2.16) in the form 



^ (Qik — qjk)Uij, M -A- + J/xjUjj, pt-r-t — U/j, M -r-lJj 



"J 



+ 0*-r-l)U«..M-r-i +*£ 2 2 B^U sj , M - r -2-A = 0. (2.21) 
For jLt=r-h2, (2.21) becomes 

2 fo* ~ ©*)Uy,r + 2-* +■ J^Uyi - Uyi J^j + U -i + B 00 = O. (2.22) 



If qi^ — qj^^^O (/a = 0, 1, . . . , r), and U#, r +2 is an arbitrary //,$ X fjij matrix, then 

(l + ftr+i-ft.r+OUtn+JMiU^-U^J^ + B^^O (2.23) 

where, under the assumptions of Theorem 2.1, <ft, r +i ~Qj,r+i ^ 1 if i ¥*j. In either case, whether 
i=j, or i ?*j, eq (2.23) is of the same form as (2.10) and hence it uniquely determines Uyi. On 
the other hand, suppose that p is the smallest integer in ^ p ^ r for which q\ v — q jf , ^ 0. Then 
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from the above analysis we know that Ugo = Uyi = . . . — \Jij, r +i-p=0 (fiiXfXj). Thus if 
Uij, r +3-p • • • , Uij,r+2 are arbitrary fiiX jxj matrices, then 

(Qip — QjpWij, r +2 - P + B ijp = O (2.24) 

and this equation uniquely determines the hitherto unspecified \Jij, r +2-p- Thus we see that the 
/i/XjLij matrices Uyi are uniquely determined for i, jf=l, 2, . . . , /, and hence that Ui has been 
completely and uniquely determined. 

Assume that for some integer ^0 we have completed solving (2.16) for /x = r+2, . . . , 
r-f2-fA; that we have uniquely determined Ui, . . . , Ui+a: and that we have determined all 
blocks of the other matrices Ui+ S , 5 > Zc, as far as possible. More precisely, we make the following 
assumptions for each (i, j)ih block Uy, i+m of Ui+ M , (0 ^ /x ^ r+2 + k). If qi^ — qj^^^O, /x = 0, 
1, . . . , r, then assume that Uyi, . . . , Uy, r +2+* are arbitrary /x*X/xj matrices. On the other 
hand, if p is the smallest nonnegative integer in ^ p ^ r for which qip — qjp ¥" then assume that 
Uiji, . . . , Uij, r +2+k-p have been uniquely determined while Uy-, r +3+fc-p, • • • , Uy, r+2 are 
arbitrary fJbiXfXj matrices, to be determined later. 



With xx=r + 3 + /c, (2.21) becomes 



k+1 / 



r+ 1 _ - { 

V (g to -^U«,r« + fc-. + Ji 4 U ,2 + *--U« > 2 + iJ^+(* + 2)U« f 2 + *+2 2 BiwU^fc+1-.^O. (2.25) 



If qis—qjs^ 0(5 = 0, 1, . , r), then with the assumption that Uy, r +3+fc is an arbitrary /XjX/Xj matrix, 
the above induction hypothesis, and the assumptions of Theorem 2.1, we see that (2.25) is an 
equation in Uy, fe + 2 of exactly the same type as (2.10) was in TV, and thus uniquely determines 
Uy, fc+2. On the other hand, suppose that p is the smallest nonnegative integer in ^ p ^ r for 
which qtp — qjp^O. Then assuming Uy, r +4+k- P , . . . , Uy, r + 3 + fe to be arbitrary /x*X /Xj matrices, 
we see that eq (2.25) may be written 



taip-yUy, r+ 3+H = ~ 2 fa* — ^)Uy,r+3+fc-«— (A+2)Utf,*+2 — 2 2 Bfc*U„j, *+i-* (2.26) 



By the induction hypothesis, all matrices on the right of (2.26) are known; hence Uij, r +z+k-p is 
uniquely determined. 

Thus the matrices U M (xx=l,2, . . .) are uniquely determined by eq (2.16). 

The above proof can be easily modified for purposes of computing the 7th column of blocks 
only. Suppose we are interested in knowing the^ coefficients of a particular vector solution, say 
the vl\\ column W*-(z) in the 7th column of blocks W/z) of W(z). If J^ . has a one in its i^th column as 
well as in its (v+ l)th, . . . , (/x — l)th but not in its /xth column where xt t=k xij, then it is readily shown 
that 

W%z)=± { I U 5 , Q 0^] *-*«*«, (2.27) 



where U s jk is the 5th (s= 1,2, . . ., /Xj) column vector in theyth column of blocks of the matrix U&. 
Thus to obtain the i>th formal solution vector in the yth column of blocks it is only necessary to 
compute the xxth, (/x — l)th, . . ., z^th columns of Uf* at each stage. 
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3. Error Bounds for Formal Partial Sum Approximations to Actual Solutions 

We shall obtain an actual solution vector of the system (2.3) to (2.5) corresponding to j>th formal 
solution vector in the 7th column of blocks of formal solutions (eq (2.27)), and we shall obtain error 
bounds for the difference between a partial sum of (2.27) and the corresponding actual solution 
vector. 

3.1. The differential Equation for an Approximation 

Starting with the partial sum 

m— 1 

*im(2)=2 U jk z- k e^ (31) 

fc=0 ' 

where in the notation of Theorem 2.2, U jk is the jth column of blocks of the matrix U* and 

€lj(z) = qjiz)!^ J Mj log z, (3 #2 ) 

we define Rj,„(z) by 

j- z ®Mz) -tO'tt + z-*X(zWjm(z) = R jm (z)e c ^ (3.3) 



On expanding the (£, ;)th block of (3.3) using (2.14) and (2.2S) and denoting q t (z) - qj(z) = qy{z\ we 
obtain the following expression for R,>,(z) 

q' u (z)V W nZ- '" - Wijiz) - (q ip - qj p )z r -P |U ijn,z~ m 

min ( M , r+l) 



+ 2 2 (qis-qjs)\Jij, ^. f + J Mf Uy, ^_ r -i - Uy, m-/-J Mj 

fi-r-2 ~\ m-i 

+ (M-r-l)Uy (M - r _,+ 2 {B,Ui ( M -r-2-,}i af-'-j {S5»-.WU>}ir^-« (3.4) 

if i 7^ ;, and for some p in ^ p ^ r, g ip - £/ p 7^ 0. On the other hand, if q ip = q jp , ^ p ^ r, then 

m-l 
R U M) = [(q>, r+l- qj, r + l + m)U Um + J^Uijn-lJutnJ^Z-™- 1 ~ ]g {85»-^)U/,} /Z"'^ 2 , (3.5) 

s=0 

where (3.5) includes the case i=j. 

3.2. The Various Possible Cases 

It is demonstrated in [4] that for each fixed j (/= 1,2,...,/) the complete vicinity of infinity 
can be divided into a finite number TV, say, of sectors *$%■ (A=l, 2, , . . , N) such that for all i 
(i = 1, 2, . . . , /) for which q {p — qj t , # for some integer p in ^ p ^ r, we have either Re qtj{z) 
bounded above, or else Re qtj{z) tending to +°° as z— ■> °° in a particular sector «5%\ This fact 
enables us to classify the various ways of constructing actual solution vectors corresponding to 
formal solution vectors under two different headings: (a) those cases for which it is possible to 
define an actual solution vector by a single Volterra integral equation; and (b) those cases for 
which we require a simultaneous pair of Volterra integral equations to define a solution vector. 

Let us establish our notation for purposes of characterizing cases (a) and (b) above. The 
fixed limits of integration in the closure, «©", of the domain 2 in which @(z) is holomorphic will 
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be denoted by {£}. In case (a) {£} consists of a single point £, say, and in case (b) {£} consists of 
two points, £1, and £ 2 , say. $ \ will denote a closed disc with centre at the origin and radius a. 
If for all integers i in the range 1 ^ i ; Si / q ip = q jp for p = 0, 1, . . . , r, then a = provided that 
every member of {£} is unbounded. , Put x = niin |£|. If x < °° we set a = x- BT» m case (b)» a ^ 

members of {£} are bounded then we assume that £1 ¥- £ 2 and |£i| = |£2|. If, for at least one i in the 
range 1 ^ i ; £ /, we have g* p ^ <£> for some p E (0, 1, . . . , r), we let a denote the maximum 
modulus of the set of zeros of q'^{z) for all such i, and set a = a if x = °°- If X < °° chose {£} in such 
a way that x = oo and no member of {£} coincides with a zero of qlj{z). Let j^* = «^ — ^ , a — ^ where 
y is empty or {oo} according to whether or not Q(z) = 0. The paths SP are defined as in [1], 

Case (a) above is related to the extreme eigenvalue case in [1]. It is characterized as follows. 
We define a region i^(z, £) to be the union of all points z, other than £ , such that there is a path 9* 
connecting zGS(z, £) and £ satisfying the following extreme eigenvalue conditions: 

(1) If £Ei^*-^*, then 9 lies entirely in ^* except for £. 

(2) For given £ and z, all elements of the matrices expjG^T) — Gi(t) — I M .[^j(r) — ^(t)]} 
(i=l, 2, . . . , /: i 7^ j) are bounded for all points t, r in the order £, t, r, z on «^. If J M . t^O, or if for 
some i(i = 1,2, . . . , /; i ^ j) we have q ip = q jp for p = 0,l, . . ., r, we require that r/t be bounded. 

(3) ^(t 17-1 ) is bounded, where the fixed number r\ may be taken to be zero if J M .=0 but other- 
wise J] is an arbitrary positive number less than 1. 

Case (b) above is related to the interior eigenvalue case in [1]. It is characterized as follows. 
We define a region i^(£i, £2, 17) to be the union of the points z for which there is a path 9 connecting 
£1, z and £2 (in that order) satisfying the following interior eigenvalue conditions: 

(1) Except possibly for £1 and £2, & lies entirely in 2*. 

(2) The integers i= 1, 2, . . . , / are subdivided into two nonempty disjoint classes Ni and N 2 
such that if £1, and £ 2 are given, then for all points t and r in the order £1, t, r, £2 on ^, the elements 
of the matrix C^(r, t) = exp {€k(r) — £k(t) — I^foM — <?j(t)]} are bounded if i E Ni (i ^ j) /while 
if i E N 2 (i 7^ j) then all elements of the matrix Qy(t, r) are bounded. If J M . t^O, or if for some 
i in the range (1, /) (i # j) we have qi P — qj P for p = 0, 1, . . ., r, then we assume without loss of 
generality that both i and j are in N 2 , that £1 is bounded, and that t/r remains bounded for all points 
t, r in the order £1, t, r, £2 on 9. 

(3) I { ai + a 4 + V (ai — a 4 ) 2 + 4a2a 3 } %> (t T,_1 )<l, where rj may be taken zero if J Mj =0,but 
otherwise 77 is an arbitrary positive number less than 1, and ai, a2, a 3 and a 4 are defined by (3.31). 

We observe that it is no longer always possible to choose the points {£} at infinity, in contrast 
to the cases in [1]. 

If r ^ we can again represent the eigenvalues of Co as points in the complex plane, and en- 
close them by the smallest strictly convex closed polygon. As opposed to the distinct eigenvalue 
case it is now possible to have several formal solution vectors corresponding to a particular eigen- 
value qoj of Co. Nevertheless, if qoj is an interior point of the polygon and if at least one endpoint 
of integration is taken at infinity 5 it will be necessary to express the error vector by a simultaneous 
pair of Volterra vector integral equations. If qoj is an extreme point of the polygon and only one 
formal vector solution corresponds to q j, it is again possible, as in the distinct eigenvalue case, 
to express the error vector by a single Volterra vector integral equation. If qoj is an extreme 
point of the polygon and more than one formal vector solution corresponds to qoj, or if qoj is any 
other boundary point of the polygon, we can always determine whether we can express the error 
vector by a single or a simultaneous pair of Volterra integral equations by testing whether or not 
the above extreme eigenvalue conditions can be satisfied. 



5 We are here attempting to construct actual solutions for which the formal solutions are asymptotic expansions. 
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We are now prepared to set up the integral equations. Let a vector Wj(z) of holomorphic 
functions satisfy the differential equation (2.3). Denoting the v\h column vector of <& jm by <£j m , we 
find from (3.3) and (2.27) that the error vector 



€] m (z) = (WjT(z) - QMtfe-W 



satisfies the differential equation 



j- z eJ m ( 2 )-[DJ(z) + 2 -^tz)]€j m ( 2 )=- £ R] m (z) 2^"" 

where Rj m (z) is the sth column vector of Rj m (z), and Di(z) = Q(z) — lqj(z). 
Hence for the extreme eigenvalue case, if the vector e"p n (z) satisfies 

where Rj*,(z) is defined by 



R"*„(z) = ~f Z e^-^> £ Rj, 



(log ty -» 

,W T7 (it, 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



and the path 0* of integration satisfies the above extreme eigenvalue conditions, then €j,„(z) simul- 
taneously satisfies (3.7). 

For the interior eigenvalue case we partition €'j m (z) into €j„ n (z) and cy m2 (z), D(z) into Dji(z) and 
D j2 (z) and R] m into Rj wl (z), Rj,„ 2 (z), and in a corresponding manner s ^3(z) into tti(z), a 2 (z), a :i {z) and 
a 4 (z) so that (3.7) may be written 



d_ 
dz 



r€ V - Az\ 



*jm2(z) 



<D' jt (z)Os i n,(2) a-,(z)\ 

WW 



\OD'(z) 



K a-Az) ajiz)/-! 



,6$ ml (z)v /RW^ ( l g z) . -,. 



(5-^)! 



v «3U f ) / 



^RJJW' 



(3.10) 



In (3.10) €/„„ (z) contains all the ith elements of ej m (z) for which i £ Ni, and €j m2 (z) contains all the 
ith elements for which i£N 2 . It follows that if e) m {z) satisfies 






Vo 1 



'e D J^)-Dj,(,) O 



'ait) ait)\/& m i(t) 



q e D#>-Ddt)/ \a : lt) ait))/\ei m Jt) 



dt+\ 






(3.11) 



^here 






(3.12) 



and the path of integration satisfies the above interior eigenvalue conditions, then €j )n (z) also 
satisfies (3.7). 
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3.3. Treatment of Ry m (z) 

We shall derive expressions for the elements of the matrices (3.9) and (3.12) which will enable 
us to obtain good bounds in the following sections. Let us consider the block 6 



RymW = - / * e^-^RijJty^-^dt. 



(3.13) 



We split (3.13) into its dominant part plus a remainder; thus we express R ijm (t) in the form 



Rijm(t)~ 



OjWUym+T U<, m -U« m Oj(t) 



r m +PisM~ 2 Wm-MVj^a-™- 2 (3.14) 



where 

PuJt) ■■ 



-[[(qt.P+l — 9j,P+l)f~ P ^ I + - • ■+(<]i,r+l — qj,r+l)t~ 1 + mt l + Jn t t l ]V ijm — UijmS^t, { ]t 



-2 



min (r-t-1, ix) 
k=/x+l-m 



+ ( M -r-l)f/ -,M-r-i+ |X 5 2 {^f/. 

A-=0 



j, M-r-2-A:} i 



t r -<" (3.15) 



The form of the first term in (3.14) is such that when substituted into (3.13) it can be readily 
integrated by parts. We define Sy(z, t) by 



Stf*,*)** 60 "*® =e Ci{z) - Ci{t) P ijm (t)e^-^ 



(3.16) 



Clearly, knowing the form of 0(0, we can express the matrix Sy(z, t) explicitly. In any case 
Sf/z, t) is simply FyM multiplied by a matrix polynomial in log (z/t). On substituting (3.14) into 
(3.13), integrating by parts and making use of (3.16), we find that 



Rft. (Z) = USJW & ~" ~ ? l ~" 1 ) + V bmW (^" m " 1 - £ T,1 " m " 1 )-|" V?. m (z, (m + 1 " ,0r<"«^l>*, 



(3.17 



where r)\ may be taken zero if J M/ — O for alli= 1, 2, . . . , /, otherwise it is an arbitrary positive 
number less than 1. The matrices \Jf jm (z), VJ- w (z) and Vf- m (z, t) are defined by 



UUz) = 



U iJm z-»'e ggg? - e c<2) " C<4) U jJm C-"'e Cja) " ?j(2) 

^r?! -m_ rr) r m 



Z^i-m-l — £fj- 



(3.18) 



(3.19) 



6 Here and henceforth we interpret e *-/') ©W = c QM Iu.9.; W- 



196 



vyz,o=- 



r"i 



m + l—ff: 



dA q' u (t) 



">2 {iB m -^)U,*}*? 



^ j (0-9j(2) 1 



(3.20) 



Equations (3.19) and (3.20) hold, provided that there is an integer p in ^ p ^ r such that qi P — qj P 
^ 0. If this is not so we define \} jm (z) = 0, and 



^7j//A 2 ? v 



*-*i 



771 ~h 1 — T}\ 



£ {Sm-^OU^iC 



£.it)-qjz) 



(3.21) 



The representation (3.17) was motivated by the following consideration. We observe that 



each element of R*- W U) is of the form 



P (z, Q=M-M)+g(z)-g(Q+ J* Hz, t)r™-*dt. 



(3.22) 



where [f(z)\=0(\z\- m \\ogz\^) (see eq (2.27)),, |g(z)|=0(|z|- m - 1 |log2^- ,/ ) and Ai(z, t) = 0(|logz| A " 
|log t\*\ as \z\ , |/|— » », where k and 5 are integers. The reason for introducing terms like t' 71 , t'^x 
where r\ and r\\ are nonnegative numbers less than 1 is to eliminate logarithmic terms conveniently, 
and at the same time obtain good bounds. From (3.22), we have 



\p(z, Q\ 



m-m 



+ 



gJz)-g(Q 






fa 7 *!-™- 1 — pV 



l ) + h(z, t)r^t^~ m ~ 2 dt 
h 



r)i -m) I 



( (r)i-m-l)(g(z)-g(0) 



V*i' 



x dt\ 



I 



{ z Vi-m-l— £ij 



+ A(z, OHi ^ 



.,}, 



2 ^ 



from which 

. sup 



\p(z, 01 



A*)-M) 



m- 



-pi- 



jr (^i-^)+ 



sup 






r (^i-™- 1 ) ( 3 - 23 ) 



where in the last supremum £ and t are points on 0* in order £, £, r, 2. 

In each of the next two sections we first establish the existence and uniqueness of the solutions 
to the integral equations (3.8) and (3.11), then we obtain the error bounds. Our procedure is, in 
one aspect, similar to that in [6] in that we first establish that the conditions in sections 3.2 deter- 
mine a domain. We then appeal to the contraction mapping principle (see e.g. [8], sec. 11) to 
establish existence and uniqueness. Proofs by use of successive approximations are given in [7]. 

3.4. Error Bounds for the Extreme Eigenvalue Case 

We use the following notation to refer to the elements of a partitioned n X n matrix A. 
\j= {A}j will denote the/th column of blocks in A, A#= {A/}f= {A}y the (i,j)th block, Ay= {Aj}" 
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the i/th vector in the jth column of blocks, and Aff= {A j }f 1J = {Afi? the (o),v)l\\ element in the (i j)th 
block. 

LEMMA 3.1: If under the extreme eigenvalue conditions of section 3.2, 2){z, Q contains at least 
one finite point z*, distinct from £, then 3)(z, Q is a domain, that is, an open connected set of points. 

PROOF: The region 3)* defined in Section 3.2 is certainly a domain. By assumption there is 
a path & in i^* connecting z* with £ satisfying the three extreme eigenvalue conditions. Simi- 
larly there is such a path connecting an arbitrary finite point z on @* with £. We can clearly choose 
a sufficiently small positive number p such that the disk ^(p): \z — zo|< p lies entirely in ^*, and 
when an arbitrary point z x of the disk is joined to z by a straight line 5?\ then conditions (1) and (3) 
are satisfied on SP U 8P\. It remains to be shown that condition (2) is also satisfied. 

To show this we define matrices 
YL&>)= sup |exp [Oi(t)-Oi(0-I|4(«(t)-«(0)]| (*»l,2 f . . ., /; i*j) (3.24) 



K/(^i) = SUP lexptCKrO-C^O-lM/^rO-^Ui))]!, (3.25) 

Ti, fi6#i 



where r, /, Ti, £i are points on ^U«^i in the order z\, Ti, fr, zo, t, t, £, and the suprema apply to 
individual elements of the matrices, as in [1]. Clearly, by hypothesis each of these matrices has 
nonnegative bounded elements. We now choose r and t on «^U^i in the order z\, r, t, £. We 
then consider 

Ki(^U^)= S " P |exp [£Wt)-0«-1i4(*W-»W)]| 

t, ^E^U«^i 

where £ is a part of the range (1, /) other thany. If r is on 3P then so is t, and in this case by hypoth- 
esis the matrix on the right consists of nonnegative bounded elements. Similarly when r and t 
are on &\. Suppose then, that r is on 3P\ and t is on 9*. In this case 

lexp [Oi(t)-Oi(0-I|4<«(t)--«(*))]H|mp [C^-CKO-Im^jUo)-^))] exp [Q,(t)-0«(*) 



- l H (qj(T) - qj(zo))] | ^ | exp [&*(*>) " &i« - I^(zo) - ©«)] 1 1 exp [Os(r) - Orf*)- l M .(^(r) - ^(z ))]| , 

(3.26) 

so that Ki(&U&i)£Ki(&>)Kl(&>i). This shows that the elements of £*(^U^i) (i = l, 2, . . ., /; 
j #y) are bounded, and therefore that 3(z> Q is a domain. 

We next establish the existence and uniqueness of the solution of the integral equation (3.8). 
To this end, we observe that all known functions in (3.8) are holomorphic and single-valued in 
every simply-connected subset of 2)(z,, £). If 3){z, Q is multiply-connected, we interpret 2{z, Q 
to be a Riemann surface on which all known functions of equation (3.8) are holomorphic and 
single-valued. 

Consider the set $${Q of all n X 1 vectors h(z) of functions holomorphic in £&(z 9 Q such that 

h(z)=0, where z^> £ along some path & belonging to the family described above. The set 

Z > b 
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,<p(£) together with the metric vector 



d(h u h 2 ) = ^|h,(t)-h 2 (t)| (3.27) 



forms a complete 7 metric space. 
Let 



US- SlWI (3-28) 



T 



e^> 



r > = *, re^ |V ^ (T)+V ^ (T ' f)l (3 - 29) 



8 = ^ 



sup 
t, tEl@ 



e W-W>«(t)r«|| (3.30) 



1-7, 

where for example Vj^ (&= 1,2) is the i4h column vector in the column of blocks \j) tl (k= 1,2), and t 
and t are on & in the order £, £, r, z. By the analysis of the previous section, all elements of the 
vectors and matrices defined in (3.28) to (3.30) are bounded. It follows that with €j m (t) replaced by 
h(t), the right of (3.8) is a contraction mapping of a complete metric space into itself, so that the 
solution to (3.8) exists and is unique. 

On substituting the bounds (3.28) to (3.30) into (3.8) and using Lemma 2 of the appendix in 
[1], we obtain: 

THEOREM 3.1: //, corresponding to a formal vector solution W[ m (z) of equation (2.3), we can 
determine a region 2)(z, Q satisfying the extreme eigenvalue conditions of Section 3.2, then equation 
(2.3) possesses an actual solution vector of functions holomorphic in 3(z, £), given by 



WUz)= 



f U il £^T ) Z " R + « <*>] f* . (3-3D 



in which 



c jm 



(z) | ^ exp {B^ (t^ 1 )} {V$n (t*«-) + y\J^ (t^--- 1 )} (3.32) 



for all z in @(z, £), where the matrix B and the elements of the vectors Uj£ a/irf y? are defined by 
equations (3.28) to (3.30). 77ie solution W[ m (z) depends on g anJ an arbitrary positive integer m. 
77ie function qj(z) is defined in Theorem 2.1. 77*e integers pi and v and the vectors U* k (k = 0, 
1,2, . . . ;s = l,2, . . . , p) are defined as for equation (2. .27). If each J M =0 (i = l, 2, ...,/) 
j/*erc */ie numbers 7) and t?i ca/i 6e toA;erc to oe zero; i/ J Mj =0 ^e/i 77 can be taken to be zero; other- 
wise 7] andrjx are arbitrary positive numbers less than 1. 

Let us now obtain a norm bound for the vector e] m (z) in Theorem 3.1. To this end we use the 
definitions (3.18) to (3.21), and define 

u SUp , , Dj(t) - Dj(t) en / x 11 

* = l^M,TE^ l|e ' ®M<-"II> (3-33) 

YS-^llujzrwn, Tj m = f? s T u e P ^ l|v ^ (T)+v ^ 2(T ' *>!!• (3.34) 



7 This is easily verified by use of Vitali's theorem. 
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where t and r are on & in the order £, ft r, 2. Substituting these quantities into (3.8) and using [1] 
Lemma 2, appendix, we obtain 

THEOREM 3.2: A norm bound for the vector €f m (z) in Theorem 3.1 is given by 

||€* m (z)p exp{S^ (t"-')}{f/^ (t"»- m ) + yr m ^ (Pi-"" 1 )}, (3.35) 

where 7/ and r)\ are defined as in Theorem 3.1 and B, U?* and y% m are defined in (3.33) and (3.34). 



3.5. Error Bounds for the Interior Eigenvalue Case 

We shall first establish the existence and uniqueness of the solution of (3.11) by use of norms. 
We split the vectors Uj^(z), \ v jm(k= 1, 2) denned in (3.18) in the manner corresponding to the way 
(3.7) was split in (3.10) into lower dimensional vectors Uj^(z), \*jm S ( s= 1»2). We change £ to £ s if 
£EN S (5=1,2: in the notation of the interior eigenvalue conditions). We then define 

fe*(!i, t>)=^ p j|vj^(t 2 )| \r u , rff 1— 1 ) 



* 2 e^ 



+ 



sup 



ft, r 2 e^ {||Vj!*(&)+VjSkte. faJlDn.^i-" 1 " 1 ) (A = 1,2) (3.36) 



where ft and ft are on 0* in the order £i, ft, t 2 , £2 if 4 = 1, or £1, ft, ft, £2 if k = 2. The fixed number 
771 may be taken to be zero if J M . = (£=1,2, . . ., /) but 771 is otherwise an arbitrary positive num- 
berless than 1. 



Similarly, we define 



a k 



sup 
ft, ftG^ 5 



1-11 



e D,,(, 2) -D,,(,.),a fc (0 



s=l, fc=l, 2 
5=2, A: = 3, 4 



(3.37) 



where ft and ft are again on 3* in the order £1, ft, ft, £2 if s = 1, A= 1, 2 and in the order £1, ft, ft, £2 
if 5 = 2, A = 3, 4. The number 77 may be taken to be zero if J^ =0 but otherwise 77 is an arbitrary 
positive number less than 1. 

Lemma 3.2: If under the interior eigenvalue conditions of section 3.2, ^(£1, £2, tj) contains at 
least one finite point z* distinct from £1 and £ 2 , then «^(£i, £2, 17) £5 a domain. 

PROOF: By assumption there is a path 3 in ^* (3* defined in sec. 3.2) connecting £1, z* and 
£2 (in that order) and satisfying the three interior eigenvalue conditions. The same path SP con- 




FlGURE 1 
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nects an arbitrary point z on 0* in this manner. By an argument similar to that used in the proof 
of Lemma 3.1, we can clearly choose a sufficiently small number p > such that when we alter 9 to 
pass through an arbitrary point Z\ in the interior of the disk S?(p) instead then (1) and (3) of the 
interior eigenvalue conditions remain satisfied. To prove that condition (2) also remains satisfied 
we need merely show that if p —» and the altered path 3P' in figure 1 tends to «^\ then 

% ^ 9 \** ' °M ++ ^r^ll^-^.^-"! (3.38) 

where t and r are points located on the paths &> and &>' in the order &, t, r, £2. A convenient way 
of achieving this is to put 

tit, t u r, n)=e l ^^- l w^aifc)rr*-« l ^ w ~ ,%,w aift)r' (3.39) 

where r and t are located as described above; outside 5%o), t\ = t, T\=t while in the interior of 
Sf(p) h and Ti are arbitrarily located on 9*\ Let us denote the norm on the left of (3.38) by a' , 
that on the right by a. Then it follows that 

a' = a+sup||A(f, t u t, t,)||, (3.40) 

where the supremum is taken for ti 9 n E ^' and MG^. The norm on the right is clearly a con- 
tinuous function of its elements, and tends to zero as p—»0. A similar argument applies to Ct 2 , 
a 3 and a 4 . Moreover the difference f^it r, ~ l )—f&( ri ~ l ) can be made as small as we please by taking 
p sufficiently small. The truth of Lemma 3.2 is now evident. 

We again observe that if i^(£i, £2, v) is simply connected then all known functions appearing 
in equation (3.11) are single valued in this domain. //, on the other hand, i^(£i, £2, r)) is multiply 
connected then we interpret «^(£i, £2, 17) as a Riemann surface on which all known functions appear- 
ing in eq. (3.11) are single valued. 

Consider the set 8 £(&, &) = (^i(£i, £2), &2(£i, £ 2 )) r of n X 1 column vectors h(z) = (h,(z), h 2 (z)) r 
of functions holomorphic in ^(£1, £2, 17), sueh that hi(z) 6 &i(£i, £2) is icXl, h 2 (z) E £) 2 (£i, £2) is 
(r-k)XI, and all elements of h(z) are bounded as z— » £1, z— » £2 along any path ^ of the family de- 
scribed above. If we take 

d(x,y)= (3.41) 

( SUP ,H X 



W)-7^)l) 



as metric for £>, where x*(f), y*(t)6Mi» £2) (*=1, 2), then it follows that £(£i, £2) together with 
the metric (3.41)iform a complete metric space. Substituting (3.36) and (3.37) in (3.11) and using 

the fact that the eigenvalue of largest modulus of the matrix I 3 )?£ (t v ' 1 ) is \ {a l J ra 4 

\a,2 •> a>4/ 

+ V(ai — a 4 ) 2 + 4a 2 a 3 }^$ (^ _1 ) and therefore less than 1, we see that the right of (3.11) is a contrac- 
tion mapping of a complete metric space into itself. Hence the solution of (3.11) exists and is 
unique. 

Following [1] we shall again first obtain a vector bound. We thus define \fj m k(u, v) and a^ as 
in (3.36) and (3.37) respectively, except that the norms in these equations are now replaced by 
absolute values. Thus \\f m k and slu are now matrices of absolute values. On substituting these 
bounds into (3.11) and proceeding as in [1] we obtain 



8 The index T denotes matrix transposition. 
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THEOREM 3.3: //, corresponding to the vector Wj m (z) (cf, eq. (2.26)) we can determine a region 
@(£ii &> v) satisfying the interior eigenvalue conditions of section 3.2, then the eq. (2.3) possesses 
an actual solution vector of functions holomorphic in <®(£i, £2, t?) and depending on m, £1 and £ 2 , 
given by 

Iwu^v Lt'oltiVu^(8-v)! r + w m2 w/J (3 - 42) 

//ere 

|€[ ml (z)| ^ exp{a 2 ^ 1 , z (t^- 1 )}i// ml (£ 1 , z) + Ci«i, zJtexpfa^'^" 1 )}^^, &) + B] (3.43) 



and ai, a 2 , a 3 , and a 4 are defined by (3.37), ^ mk (u, v) (k= 1, 2) are defined by (3.36) with norms in 
these equations replaced by absolute values, and 

B= [l 2 -C 2 C^ u fo)]^C 2 X exp {aft (t^)}V ml (tu «' + C,({i, &) exp {a 4 ^ (t^ 1 )}^^!, £2) 

(3.44) 



Ci(fc, z) = a 1 - 1 [exp {a 1 n ] ,^- 1 )}-I 1 ]a 2 ; C 2 = a- 1 [exp {a 4 ^ (t^ 1 )} -I 2 ]. (3.45) 

In (3.44) Ii designates a k X k unit matrix and I 2 designates an (n — k)X (n — k) unit matrix, k being 
the number of elements in Ni (sec. 3.2). All variations on the right of (3.43), (3.44), and (3.45) are 
taken along 8P. The function qj(z)< is defined as in Theorem 2.1. The vector U s given by Theo- 
rem 3.1 is partitioned into lower dimensional vectors as described in section 3.2. The numbers 
7] and 171 are' defined as in Theorem 3.1. The bound (3.43) is valid when every eigenvalue of the 
matrix C 2 Ci(£i, £2) is less than 1 in magnitude. An exactly similar result holds for € v . m2 {z). We 
also have 

Theorem 3.4: Let ai, a 2 , a 3 , and a 4 be defined by (3.37), ^ mk (k= 1, 2) be defined by (3.36), and 

(3.46) 
F ( z ) = (aXlif [e(ai+H4)f(Z) " l ~ (ai + **&*)] ' 



where the variations are taken along 3P. If F(£ 2 ) < 1 then 1 — F(£ 2 ) times the norm of the vector 
€ jmi( z ) ( c f- ea - (3-42)) is bounded by the expression 

e M«[F(z)+l-F(&)]^ m i(Ci,z) 

+ a2&3 [l- e -< a . +a ^ z »] [e< a * +a <> f <y-e< a > +a <> f < z >]^ ml (z, &) 
(ai + a 4 ) 

+ -^— e a .« z » + V(W[l -e-^ + ^ f ^]^ m2 (Cu &) (3.47) 
ai ~r a 4 

An exactly similar result holds for € P im2 (z). 
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4. Example: The Case n = 2 

In this section we consider the special system X' = z h ^l(z)X (h an integer), that is 

An *J*\ (a u (z) a V2 (z)\/xn x i2 \ 

Y*2i x 22 J \(Hi(z) a 22 (z) \x 2l x 22 J 



where for each arbitrary positive integer s 

s-l 

«oU) = 2 a tikz~ k + ao.vU)z- s (*, j= 1, 2), (4.2) 

A*=() 

and the functions ajj. s (z) (5 = 0, 1, 2, . . .) are holomorphic in a domain Sf 1 and uniformly bounded 
on some path & E 2' extending to infinity. 

We begin by transforming the system (4.1) to canonical form. The flow chart on the following 
pages indicates the steps involved in transforming the system (4.1) to canonical form. The chart 
is constructed along the lines of the proof of Turrittin's theorem [3], and it contains all the trans- 
formations used in [3] for arbitrary n. The process terminates in all cases. 

Box #2 in the chart contains the normalizing transformation, box #8 the exponential trans- 
formation, box #19 the shearing transformation, box #27 the zero-inducing transformation, and 
box #31 the root-equalizing transformation. Box #7 also contains a modified form of the zero- 
inducing transformation. 

The resulting system takes the form 

r/W 

— = z'U~(z)W = [Q'(z) + z" 2 «(z)]W ( 4.3) 

in which r is an integer, and 



*v*4 m ° W° %A( bM M 1 

l\0 qfc)/ \f*-> 0/J lW(r) b 22 (z)) 

0'(z) is the first matrix in braces on the right, 93 (z) is the second, and 

qtk 



(4.4) 



*(*>= 2 JV-t * r+ '~* + g'."-' h % z (i = L2). (4.5) 



fc=0 



If r ^ — 2 then C(z) = 0. If for some integer p in the range ^ /; ^ r+ 1, we have </i ; , # ^, then 
£ = 0. Otherwise f is either or 1. The functions bij(z)(i= 1,2) have the expansion 

bij(z)= J 6^**-*+ bi jt (z)z-; (4.6) 

A*=0 

where the elements of 6yt(z) are holomorphic and uniformly bounded on a path ^ extending to 
infinity in 2. 

By the analysis of sections 2 to 4, the system (4.3) has a formal solution matrix of the form 

W«=U«e°w (4.7) 

203 



FLOW CHART ILLUSTRATING TRANSFORMATION OF THE 
SYSTEM TO CANONICAL FORM WHEN n = 2 




Make transformation 
X= T Z taking A 
to Jordan canonical 
form. 




Yes 




Yes 



Make transformation 

X = [l + T ] z _1 + T 2 z~ 2 +- — + T r+1 i" H ]z 
diagonalizing A 1 , A 2 , • • , A f+1 



V 



■© 



Make transformation \ * 

X=Ze Xz with(r+l)X 
eigenvalue of A . 




G^D 



10 




17 t 



>— © 



Figure 2 
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21 
Make transformations 

r * 2r 

'ijk 
b ij,2k+l^ 0{i #j); 

b ii, 2k * — 0; 

b 12,2k M 2a 12, k*l 

b 21,2k * 2a 21,k 

b ii,2k + l * 2a ii,k+l 

(k= # 1 # 2, — ) 



l + a , 



u ll, 2r + l ' ' u ll,r + l 

a ijk * b ijk 

[k= 0,1 # 2,---J 



24 



Eigenvalues X^ , X2 define 

integers k 1 , k 2 , * and 

number s by 

k 1 < Re \ } < kj +1 
k 2 < Re X 2 < k 2 +1 

k~* — I Ic 1 -k 2 , s-<— X 1 - X 2 



31 



If solution corresponding to 

X 1 i s c 

motion. 




22 



Redefine a-^ 's 
(i/j) and a n ^ f + ] by 



>21k 



"21,k-fL 



b 12k ** a 12,k+/x 

(k«>*l,-M*2 f --) 

a ijk * — b ijk (i*j) 
(k = - M +l,- M +2,---) 



'll.r + l 



M + a ll,r + l 



-0 







23 




\f 30 



X 1 is desired, make transfor- 



z h+V 



X = 



Otherwise 



ko+\ 



GE°MI> 



32 



<2> 



Apply inverse of each 
exponential transfor- 
mation made. 



— ►Tend] 



Figure 2— Continued 
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where 






uw^fe"'; U « =I 

fc=o 

(Ullfc «12*\ 

U = 0, 1,2, . . • 
U21fc U22fc/ 



Suppose that p is the smallest integer in the range Ogpgr+1 such that q lp * q ip . Then 
the iMjk are computed from 9 



lino =1 (i — 1 , 2) 

uy fc = if i^j and 0^/c^r+ 1-p (i, .7 = 1,2) 



Uiifc 



=-±'2 t (biuUu,k-i-*+t>a,U2i,k-i-,) (i = l,2;ASl) 



(4.9) 



1 r+1 

Uij,r + 2 + k-p = ^ lMj,r+2+k-s(qjs-qis) + (k-\-l)lMj,k + l 

Qip QjP s =p+l 



k 



If, on the other hand, <f lp -</2p = for p = 0, 1, . . ., r+1, then U ft+1 (ft 2 0; U = D is com- 
puted from 



(ft + 1)U* + , + JU fc+1 - U fc+] J + 2 B S U*_ S = O (4.10) 

with 

J= /° °\ £= orl. (4.1D 

In either case, whether £ = or 1, eq (4.10) with Uo — I uniquely determines each U*, A: = 1. If 
£ = each column vector can be computed independently of the other. If f = 1 the second column 
vector in Ujt (A; = 0,1, . . .) can be computed without computing the first. If, on the other hand, 
f -=1 and only the first column vector in U* (A: = 0,1,2 . . .) is required, the second column vector 
must also be computed, and the computations for the second column vector must precede the com- 
putations for the first. 

We shall obtain a bound on the difference between a partial sum of the first formal solution 
vector and a corresponding actual solution vector. Since the analysis is similar to the preceding 
extreme eigenvalue case, we only present the main points. 



9 The sum V ' s to be replaced by zero if p = r + 1 . 
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In order to define a domain, we choose a nonnegative number a such that if q\p = qip for 
p = 0,l, . . ., r, then a = if |£|=°° and a=\t\ if |£|<°°. Let us denote the largest zero of q! 2l (z) 
by ao. If for some p£(0,l, . . ., r) gi p t^ q 2p then we set a=|ao|if |£|=°°, while if |£|< °° we as- 
sume that £ E 2$ , £ does not coincide with a zero of g^z) and |£|^|a |. Let ^* = ^ — ^ a — S where 
g'a is the closed disk with center at the origin and radius a, and S is empty or {<*>} according to 
whether or not €l'{z)=0. We define a region 2(z, Q to be the union of all points z, other than £, 
such that there is a path 8P connecting z and £ which satisfies the following conditions: 

(1) If £Gi^*-^*, SP lies entirely in S* except for £; 

(2) For given £ and z, exp {^2i(r) — q>\{t)} is bounded for all points t, r in the order £, £, r, z 
on ^. If q\ p = qip for p = 0,l, . . ., r and either gi, r +i t^ 92, r+i or f^O we also require that r/t 
remain bounded; 

(3) Yg> (P*- 1 ) is bounded, where 77 may be taken to be zero if f = 0, otherwise 77 is an arbitrary 
positive number less than 1. 

By section 3.3 3f(z, £) is a domain. 

If q 2 i(z) = q 2 (z) — qi(z) ?* (and hence f = 0) let 



Rllm(z) = mi4um2; m 1_ ^ ['«. m-*C*) tt ll* + &12, m-*(z)Hl2jz 



R21m(z)= ^l(z)tt21mZ m — 21 (q2k-qlk)z r k U 2 \mZ 



k=p+\ 

m + r+l ["min (/x, r+1) 
— ^ 2 (<?2s — gi s )w 2 i, M - s + (/X — T— 1)^21, M -r-l 

IA=m+p+l ^-s^p.+ 1-m 

p-r-2 -| m-l 

+ "V (^lls^H, M-r-2- S + 6l2»W21, n-r-2-s) Z r ^~ ^ [ *"• m-fc)u>lU+ &12, m- S (* 2 , Jz" m " 2 , 
.±0 J *=0 



(4.12) 



where p is the smallest nonnegative integer such that q 2p 9^ q lp . On the other hand, if q 2 {z)= qi(z) 
then 



(Rilmiz)) 
\B.2lm(z)J 



'Hum Ui2m\{0 0\ lm 0\ / U\\ m U 12 n 
\U 2 im U22m) \£ 0/ \£ ml \U2im U 2 2m 



bn,m-s(z) bi2,m-s(z)\ fUns Uus) 
K b21,m-s(z) bl2, m-s{z) \U 2 \ S U^UJ 



\i log zj 



• (4.13) 



The integral equation for the difference between an actual solution vector and a partial sum 
of a formal solution vector thus takes the form 

€i m (z)= f*exp [€l(z)-£l(t)-l(qi(z)- qi(t))]{r 2 Wt)*Ut)-'RUt)}dt (4.14) 



where the path of integration satisfies the above conditions, and Ri m {z) = (RnmCz), R 2 imCz)) T . 
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We define nonnegative numbers uf lm and u% Xm to be the suprema when rG# of the respective 
expressions 



(ttllm + gttl2m ^g T)T~ m — (u nm + g^lm fog Qj~ m 



T i) x -m _ £17 



(«21m + ftt22« log T)T- m ^ 1 (r)-fei(i)|-| tt21|n +^ 



log - (linn, + W 2 lm log J) + U 2 2m log J 



T i7 1 -m_ £17 



(4.15) 



where rji may be taken to be zero if f = but otherwise 171 is an arbitrary positive number less than 1. 
Similarly with t\ and £ 2 points on 3P in the order J, h, t 2 , z we define 



ynm = 



sup 



*G ^ 



^-t)! m-1 

T i _ 2 {*U.»»-*(0[Hll*+f l<>g* ^12s]+ 6i2,m-sW["21 S +^l0g^ II22*]} 



721m = 



sup 



P2l(t 2 )l gglfe) -gg2itf2)-g2i(0 p 21 (g)/ g^(Q | r r? 

qi-m-Tgjj-m 



' 1-171L 



m+ . 



*1\ <?2l(*l 



*m+2 



+ 



2 ^ {U2*, »-.(*!)+ flog (~) ilM-*i)l [**U + £ logtHtej] +*P*-j-( B ¥^) 



(4.16) 



where 



m— 1 

P2lm(t) = R 2 im(t)-q 2 i(t)U2imt- m -\- ^ [2>n, m-8(t)u 1U + b 12 , m-s(t)u 2 i Jr m " 2 (4.17) 



with R 21m given by (4.12), and t u t 2 points on & in the order £, fc, fc, 2. Note that p 2 i(*) = if the 
smallest integer p such that <j 2p ^ <?ip is r+ 1. 
We define 



(l — rj)a k = 



sup 



r* b lk (t)\ 



(1 - rj)a, 



sup 



t^ 



M*i)+f log ( j) M*i)l e*i«*>%A> 



(4.18) 



(A=l,2 ; 5 = A + 2). 



In (4.18) *, and * 2 are defined as in (4.16), and 17 may be taken zero if f is zero; otherwise 77 is 
arbitrary positive number less than 1. 



an 
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Summarizing we have 

THEOREM 4.1: //, corresponding to the formal vector solution Wi(z) of equation (4.3), we can 
determine a domain @(z, Q satisfying the conditions above, then the equation (4.3) possesses an 
actual solution vector of functions holomorphic in 3)(z, £), given by 



where 



W lm (z) = 



|€llm(z) 



W llm (z) 
W 21m (z) 



^Ui2k^ 



k -°l\U2lJ W/| 



f€n m (z)\ 
V € 21m(zV 



»<?!(z) 



fan 



a 3 \ 



Urn 



'(t"'- m ) + y„ m 3£(t*i-*-i)l 



|€2i m (z)|/ [\a 2 a 4 / J [ u * m ^ (t^- m )+y 2 im^(t T,1 - m - 1 ) 



(4.19) 



(4.20) 



for all z i/i ^(z, £). 7%e solution Wi m (z) depends on £ and an arbitrary positive integer m. 77ie 
function qi(z) is defined by (4.5); */ie numbers Uj* m , yum ana 1 ai, a 2 , a 3 , a^ are defined by (4.15), (4.16), 
and (4.18). If f = 0, ooJ/i t/ and 171 may 6e £aAezi to 6e zero; otherwise rj and r\\ are arbitrary positive 
numbers less than 1. 

The right side of (4.20) may be explicitly evaluated by use of 



exp 



{;?}■ 



sinh x l~n~ fi 



8-a\ 



ih X l] 



e l/2(a+8). 



x=* 



V(a-8) 2 + 4/3y. 



5. The Choice of Paths 

The conditions on the paths of integration in this paper are somewhat weaker and define a 
larger domain of validity for the asymptotic expansions than the conditions given in [1]. The paths 
which minimize the error bounds under the present conditions are not known. In application, we 
suggest that paths similar to those of [5] be tried; the possible loss in accuracy is offset by the 
greater ease of evaluation of the variations. 

Often the error bound can be appreciably reduced by the following simple procedure. Con- 
sider for example, the distinct eigenvalue case 1. Let T be a matrix such that T^AoTo = A is 
diagonal. Then (T () D) _1 Ao(T D) =A, where D is an arbitrary nonsingular diagonal matrix. We 
choose D so that the least upper bound of the off-diagonal elements in the resulting coefficient 
matrix is minimized on the path 8P. This applies also to the more general case of this paper, ex- 
cept that here we may destroy the similarity of the canonical form to the Jordan canonical form. 
If we want to preserve the canonical form then for every Jordan leading subblock of t£(z) of the form 

r 

2 qikZ r - k l+ qi, r + J log z 

where J has a full set of units in its first lower diagonal and is zero elsewhere, we must make the 
corresponding elements of D equal. 

On inspecting the proof of Theorem 2.1, [1], we observe that the undeterminable diagonal 
elements of T*(A: = 1,2, . . ., r+1) were chosen to be zero. By introducing nonzero diagonal 
elements we may again be able to improve on the error bounds. A similar observation applies 
to the zero-inducing transformations used in [3]. 
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6. Summary 

The results of this paper and the preceding paper extend the work of Olver [5] for the second- 
order equation to a system of first-order equations. In the majority of applications it is expected 
that the eigenvalues of the lead coefficient matrix A are distinct. The analysis for this case is 
considerably simpler and different from the analysis for the case when the eigenvalues for the lead 
coefficient matrix are not distinct, and this is the main reason that the presentation has been split 
into two parts. 

In all cases the procedure was to transform (1.1) to canonical form, to obtain formal solutions 
for the transformed system, and to obtain actual solutions such that the formal solutions are 
asymptotic expansions of the actual solutions of the transformed system, as z— >°o in some sector 
of the complex plane. Two cases arise when the eigenvalues of A are distinct: the extreme eigen- 
value case and the interior eigenvalue case. They are characterized geometrically by the relative 
location of the eigenvalues of A in the complex plane. The main distinction between these two 
cases is that in the extreme eigenvalue case it is possible to express the error vector — that is, the 
difference between an actual solution vector and a partial sum of a formal solution vector — by use 
of a single Volterra vector integral equation, while in the interior eigenvalue case it is necessary to 
use a simultaneous pair of Volterra vector integral equations. As a consequence, the error bounds 
are considerably sharper in the extreme eigenvalue case. 

When the coefficient matrix in (1.1) is regular for all sufficiently large z and when the eigen- 
values of A are distinct, it is always possible to choose end-points of integration in the Volterra 
vector integral equations at infinity. This is no longer the case when the eigenvalues of A are 
not distinct. Further, in this latter case we can no longer always give a simple geometric char- 
acterization of extreme and interior eigenvalue cases in terms of the relative location of the eigen- 
values of A . We can, nevertheless, always express the error vector by at most a simultaneous 
pair of Volterra vector integral equations. 

The case n = 2 is likely to be the most important in applications. The complete solution for 
this case is given in detail, including a full statement of the transformation to canonical form. 
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